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Abstract 



We investigate a system coupling the nonlinear Schrodinger equation and the inviscid 
Burgers equation, which models interactions between short and long waves (for instance in 
fluids). Well-posedness for the associated Cauchy problem remains a difficult open problem 
and so we tackle it here via a linearization technique. We establish a linearized stability theo- 
rem for the Schrodinger-Burgers system when the reference solution is an entropy-satisfying 
shock wave to Burgers equation. Our proof is based on suitable energy estimates and on prop- 
erties of hyperbolic equations with discontinuous coefBcients. Numerical experiments support 
and expand our theoretical results. 

1 Introduction 

1.1 Background 

In [5] , Benney introduced and studied several systems of partial differential equations vifhich model 
the interaction of short waves and long waves arising in a variety of physical applications. These 
models are relevant, for instance, for the description of gravity waves in fluids or internal-surface 
waves. Several such systems were studied in the literature, especially systems that introduce a 
coupling between the nonlinear Schrodinger equation (for the short waves) and a general partial 
differential equation (for the long waves). In the literature, the latter is often a linear transport 
equation or a nonlinear hyperbolic balance law (which might also include dispersive or dissipative 
terms). More precisely, the general class of Benney's systems reads 



where Ci,C2,q;, /3,7,^ and v are real constants, and the unknown functions u — u{t,x) and v — 
v{t, x) represent the short waves and the long waves, respectively. For instance, Tsusumi and 
Hatano [25, 26 studied the coupling with a general linear equation in v, while Bekiranov, Ogawa, 
and Ponce J7i studied the coupling with the Korteweg-de Vries (KdV) equation. More recently, 
there has been a renewed interest in the (particularly challenging) coupling with nonlinear balance 
laws, which motivates us here to consider the following nonlinear Schrddinger-inviscid Burgers 
system (with e > 0): 
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i dfU + icidxU + dxxU — auv + ^\u\'^u^ 
dtv -I- C2dxV -I- ^xdlv + vdxv"^ = I3dx{\u\^), 



(1.1) 



iUt + Uxx = vu 




(1.2) 



This is clearly an important prototype for the class of systems (1.1 1. In [TF, the existence of 
local-in-time and smooth solutions to this system was established by solving the initial value 
problem, but the possibility that smooth solutions blow-up for sufficiently large times was left 
open (although some partial results are available in [5 ). On the other hand, the existence of 
global-in-time solutions to the corresponding initial value problem is known [TT' only in the class 



of weak solutions, provided the nonlinear flux v in (1.2) is replaced by the cubic function v (or 



more generally, by any concave/convex function). 

Recall also that Dias, Figueira, and Frid 13] (cf. also [UlIIS]) treated a class of Benney-type 
systems in which an additional coupling function was introduced in order to "tame" the nonlinear 
coupling between the two equations in (jOJ. These authors treated long waves modeled by, for 
instance, the Navier-Stokes equation or the p— system of nonlinear elasticity, and established a 
well-posedness theorem for the Cauchy problem. Their theory encompasses only variants of the 



Schrodinger-Burgers system but, due to the presence of the additional coupling function in (1.1 ), 



the Schrodinger-Burgers system (1.2) of interest in the present work is not included. Finally, we 
refer to 1 for the coupling with nonlinear viscoelasticity and to |31 HI [TU] for the numerical 
analysis of short waves-long waves problems. 

1.2 Objective of this paper 

Therefore, despite recent efforts on this subject, the existence of weak solutions to the Schrodinger- 
Burgers system remains an open problem. In the present work, we shed some light on the systems 



( 1.1 ) via a linearization approach and, specifically, we establish that solutions to ( 1.2 ) are linearly 
stable, when the reference solution consists of an entropy-satisfying shock wave to the inviscid 
Burgers equation. Observe that the linearization of nonlinear hyperbolic equations (and systems) 
in a class of weak solutions was first investigated in LeFloch [20l [21] . More recently, this question 
was revisited by Godlewski and Raviart jl6j from a numerical analysis standpoint and for systems 
in several space dimensions. 

Schematically, we will adopt the following approach. Given a system of partial differential 
equations and a solution U = U^{t,x) which we assume to depend upon a (small) parameter 
(5 > 0, we can formally seek for an expansion of the form 

[/^ = [/(o) +,5C/(i) + ..., (1.3) 

where U^^'^ is the chosen reference solution. Formally at least, one can linearize the given (nonlin- 
ear) system around t/*-*^-* and derive a system for the first-order perturbation [/^^^. If we manage 
to establish an estimate of the form (on each compact time interval [0,T], say) 

llC/W(i)||<GT(||C/(^'(0)||) (1.4) 

for some norm || • || and some continuous function Gt ■ ^+ IR+ vanishing at the origin, then we 
may expect that the nonlinear stability estimate \\U^{t)-U'^°'>{t)\\ < Gt(||/7''(0)-C/(°)(0)||) would 
then follow. This classical strategy was followed in LeFloch |2Ij for solutions to linear hyperbolic 
systems with discontinuous coefficients, and estimates in the norm or in the total variation 
were derived therein, leading indeed to a proof of the continuous dependence of solutions for 
nonlinear systems of conservation laws. Concerning linear hyperbolic equations with discontinuous 
coefficients, extended work is now available and, in addition to to [201 IHl US] , we refer the reader 
to [351 [II] as well as Bouchut and James [S] , who introduced a notion of duality solution for scalar 
equations. 



For the system under consideration in the present work, the nonlinear stability estimate (1.4 1 
remains a very challenging open problem so that the present work focuses on a linearized anal- 
ysis which, as we will see is already quite involved. We are going to analyze solutions to the 



Schrodinger-Burgers system ( 1.2 ), when the reference solution consists of a shock wave (in v) and 



a smooth oscillatory traveling wave (in u). Our main result will be a rigorous derivation of a 



stability estimate like (1.4). 



2 



1.3 Outline of this paper 



To start off our anafysis, we need to describe tlie particuiar solution, wliich will serve as a reference 
solution. For clarity in the presentation, we consider first the limiting case e = 0. It is not difficult 
to check that there exists a (stationary) traveling wave solution 



{u,v){x) 
(p{x) 



(e^**r(x),^(x)), 
-£.sgn{x), 



X e 



in which b is an arbitrary real parameter and (/? is a stationary shock wave (normalized, from 
now on and without loss of generality, so that ^ = 1) to Burgers equation, while the function 
r e W'^'°°(U.) is a solution to the ordinary differential equation 



— r" = — 6r + r sgn x. 

We look for solutions that are continuous with continuous derivative at x 
obvious notation, 

r{x) = ta:>or+ix) + l^<or_(a;), 
and ^+(0) = ^-(O) and r^(0) = r'_{0). Elementary methods lead to us, ii b < —1, to 



(1.5) 

0, so that, with 



r+(x) = C 



|1 + 6| 
1 - b 



sm 



(Vl - bx) + A cos {Vl - bx), 



r_{x) = C sin (^|1 + b\ x) + A cos (^|1 + b\x), 
[ r(0) ^ A, r'(0) = C^TT^, 



and, if —1 < 6 < 1, to 



(1.6) 



r+{x) = A 

r„ (x) — Ae^ 
[ r(0) = A, 



6 + 1 

1-6 

b+lx 



sm 



(Vl ~bx) + A cos (Vl ~bx), 



r'(0) = AVl + 6, 



where A, C are arbitrary constants. 

More generally, when e > 0, the corresponding solution ri,(x) is given by 

r" (x) = br^{x) — r^{x) sgn x \/ er'^{x) + 1 — er1{x), 
re(0)=r(0), r^(0)=r'(0), 



and 



(^e(x) 



sg: 



:na; Ver2(a;) + 1 G L°°(M). 



(1.7) 



(1.8) 



(1.9) 



No closed formula is now available but standard arguments for ordinary differential equations 
(i.e. Picard's theorem and an energy estimate) yield the existence and uniqueness of a global 
solution e M^^'°°(]R) of (1.8 1. An argument of continuity with respect to parameters shows 
that, as e — > 0, one has — > r in Cjqj,(M) (uniform convergence of first -order derivatives), where 
r is defined by (1.6)-(1.7). 

Obviously, the choice of a minus sign in ( |1.9[ ) is not arbitrary. Indeed, one could attempt to per- 
form the foregoing analysis with the (monotone increasing) weak solution (pe{x) = sgn x er'^ {x) + 1. 
However, as we shall see, the choice ( 1.9 ) correspond to an entropy-satisfying shock and guarantees 
uniqueness to the linearized problem. (Infinitely many solutions would otherwise be available.) 
Hence, the system (1.2 1 may be viewed as a diffusive regularization of Burgers equation which 
properly takes into account the small-scale effects, and thus allows us to consider the whole model 
as a relevant physical model. 
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An outline of the paper follows. In Section[2j we consider the weakly coupled problem obtained 
by setting the parameter e to zero. This regime is physically interesting, since (as was argued in 
[Uj) it reflects the fact that a shock wave is a macroscopic phenomenon while the Schrodinger 
equation models a microscopic one. One main difficulty stems from the fact that the linearized 
system involves a transport equation with discontinuous coefficient. In the case e = the linearized 
hyperbolic equation can be solved explicitly and yields a singular solution containing a Dirac 
measure |20) . Our main task is to analyze the interaction between this measure-solution and the 
solution to the Schrodinger equation. 

In Section[3j we consider the full problem with e > and we establish our main stability result. 
The analysis of the linearized system in this general regime is made challenging by the fact that 
the Schrodinger equation now has a nonlinear right-hand side coupled to the hyperbolic equation, 
and it is convenient to rely on a change of variables introduced earlier in |151 [1] in order to derive 
a formulation amenable to Kato's semi-group method [TH]. Moreover, the linearized hyperbolic 
equation still involves the product of a measure by a discontinuous function and can no longer be 
solved explicitly. It is more delicate to define a suitable notion of solution. Finally, in Section [4j 
we present some numerical experiments and illustrate our theoretical results. 



2 Linearized stability for a weakly coupled system 
2.1 Stability statement 



Our first objective is the linearized analysis of the system (1.2) when the coupling parameter e is 
set to zero. Thus, we consider the weakly coupled system (t > 0, a: € M) 



iUt + Uxx = vu, 

vt + {v^)x = 0, 



(2.1) 



We linearize this system around the solution (e'''*r, (p) described in the introduction, especially 
(1.6) and (1.7). We make the substitution 

V i-i' ip + Sv, M n- e'''*r + (5e'^*u 

for some small parameter S > 0. By identifying first-order terms in S, we find 



iut + Uxx — iv + b)u + vr, 
Vt + 2{(pv)x = 0, 

which we supplement with the initial data 

u{x,0) = uo{x) e H^{R), v{x,0) = vo{x) £ H^{R). 

(The notion of weak solution used in the following theorem will be clarified shortly.) 



(2.2) 



(2.3) 



Theorem 2.1 (Linearized stability for a weakly coupled system). Let r — r{x) be defined by (1.6 1 
or ( |1.7[ ) and ip — — sgn(x). Then, the solution (e*^*r, (p) to system (2.1) is linearly stable in the 



following sense: for any T > and any initial data (2.3), the system ( |2.2[ ) admits a unique weak 
solution {u,v) = {u{t,x),v{t,x)) satisfying 



where Gt 



||(u,Z;)||Loc(o^T;_ffix//-i) < GT(||('"(0),w(0))||Hix(LinL2)), 

^ is a continuous function vanishing at the origin. 



(2.4) 



The solution v in Theorem 2.1 was first introduced in [20] and can also be interpreted as the 
unique duality solution to the second equation in (2.2) in the sense of [5]. 
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2.2 Proof of stability 

We now discuss the relevant notion of weak solutions and we establish the linearized stability 



property of interest. We need to establish existence and uniqueness for the problem (2.2|-(2.3| 



and derive the bound (2.4). First of all, we observe that the second equation in (2.2) admits an 

], even for general Riemann data); 



explicit solution (determined in 



v(t,x) ^v{t,x) + '^{t)6^^, 



(2.5) 



where 



v(t, x) 



vo{x — 2t), X < 0, 
VQ{x + 2t), X > 0, 



2t 



vo{x) dx, 



and Sj2t denotes the Dirac distribution concentrated on the line segment Et = {(0, t) : t G (0, T)}. 
Thus, if is a test function in R^, we set (^Sy,0) = J^^ 4>{0,t) dt. In all that follows, and since 
the final time T > is arbitrary but fixed, we simply use E to denote Ey. 



Some remarks about ( 2.5 ) are in order. First, observe that v is clearly a solution to (2.1) outside 



the line E^. Second, observe that replacing the expression of v in (2.1 1 leads to a product of the 



discontinuous function 2^p by the distribution S^, which has no meaning in distribution theory. 
Therefore, in agreement with the definition introduced by LeFloch [50], one defines 2Lp6s = cr, 
where a is the speed of propagation of the shock ip which coincides with the direction of E (in this 
case, (7 = 0). This gives 

vt + 2{ipv)^^v + 2{ipv)^ + ^'it)6^ + ^{t){dt+ad^)6^, 

so that the definition of the product 2ipSY: merely translates the natural fact that Sy: is invariant 
in the direction of E. With this in mind, it is easy to see that Vt + 2{ipv)x = 0. 
Let us now show the estimate 



ll«IU^(R+;if-i(R)) < c(||i;o||l1(R) + |ko||L2(R)), 



(2.6) 

which will establish the linearized stability for the second equation in ( |2.l[ ). First, it is obvious 
that < ||wo||l2- Next, given (f) = 4>{t,x) e L^{0,T- H^), we can write 

\{^{t)6^A)\ - \{SY,^{mt,x))\ 

/ / VQ{x)4>{Q,t)dxdt <\\vn\\Li ||0(t)|Uoo(R)dt. 

/O J -2t Jo 

Since iJ^(M) is continously embedded in L°°(IR), we obtain 

|(*(Ofe,0)| <c||t;o|lLi|l0||Li(RH.;Hi), 



whence (2.6 1 



We now turn our attention to the Schrodinger equation in (2.2). In order to establish the 



existence of a solution, we regularize the singular term v r and, for each n > 0, we consider the 
problem 

u"(x,0) = <(x), 

where 

h"{t, x) := v{t, x) + ^{t)np{nx), (2.i 



(2.7) 



where p is a standard moUifier, V{x) = Lp + h, f = — sgnx, and v is given in (2.5). In addition, 
Uq G i7^(E) is some regularization of the data uq, but our final estimate will actually hold with 
uo e H^, only. 

Standard techniques (cf. for instance |5]) about the semi-group generated by z(A — V) with 
a potential V e L°° imply that, for each n, the problem (2.7) with Ug G admits a unique 
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solution u" e C^{[0,T];L^) D C{[0,T];H^) for all T > 0, provided that (for fixed n > 0) /i" € 
W^1'^(0,T;L2(M)). This latter property follows easily from the fact that dtK = dtV + 2{vo{2t) + 
vo{—2t))np{nx), so that (for fixed n) 

\\dthn\\L^a,T-L^) < c||i'o||l1(0,T;H1)- 



In oder to justify the passage to the limit n — > oo in the equation (2.7), we need first obtain a 
(uniform in n) estimate for /i", namely 



ll^n||L~(0,T;//-i(R)) < c(||uo||l1(r) + Hi'O || (R) ) 



(2.9) 



for some uniform c. Given </) = (j){t,x) G L^{0,T; H^), by neglecting the term involving v (which 
was already dealt with) and by using the unit integral property of the mollifier, we find 



\{-^{t)np{nx),(l)it,x))H-ixm\ < 

< IK'olU 



2t 



|tio(s)| ds np{nx)\(/){t, x)\ dxdt 



-2t 
T 



.dt<c\\v^\\Li / UWn^dt, 



from which (2.9) follows. 

Next, we derive an estimate of the L^-norm of m". Consider the equation (2.7), multiply it by 
m", and take the imaginary part: Im(wut + uuxx) = iTniurhn), where, from now on, we omit the 
superscript n. Since Ivaiuut — |i9t|up, by integrating over M, we find 



ruhn dx, 



where we have used Im {u^l'^ — 0. Integrating over (0,T) yields us 



\u{t)\'^dx <c + 









f 


/ ruhn dx 


dr < c + c / 


Jo 


Jr 


Jo 



which, by (2.9), gives 



||w"|li2 <c + c f ||M"|Uidr. 
Jo 



(2.10) 



We now estimate the gradient of u", as follows. Multiply the equation (2.7| by Ut, take the 
real part, and integrate over M in order to obtain 



1 d 
2di 



\u,\^ + V\u^\dx 



d_ 
dt 



rhn Re u dx 

rhnileudx+ I rdthnKeudx. 



Then, an integration over (0, T) gives (recalling that r is bounded) 

1 

2 



\ux\^dx+ - I Vlu-'ldx 



<c + c 



/ hn Re u dx 




/ dthjiKeudx 


dr. 


Jr 


Jo 


Jr 





(2.11) 



We now estimate the terms on the right-hand side. First, using (2.9) and a weighted Young 
inequality, one has 



/ 

Jr 



hn Re u dx 



<c{\\v,\\h + \\vo\\h) + \\\u\\H 



(2.12) 
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and, in view of (2.8), dthn = dtv + 2{vo(2t) + VQ{—2t))np{nx) , we find 



dt hn Re u dx 



dr < 



< c 



Vtu\dx + c / \vo{2t)+vo{—2t)\np{nx)\u\dxdT 



dr. 



[voWhAMl^ + \\vou\\ 



Since <Z L°°, it follows that 

ft 



dthn Reudx 



dr < 



Next, using (2.12) and (2.13) in (2.11) gives 



-^\\<\\h<c + c\\u-\\l.+c 



\vo\\hAMh^ dr. 



\u"\\hi dr, 



(2.13) 



which, along with (2.10), implies < c + c dr. An application of Gronwall's 

lemma then leads us to < c{t) for some continuous function of t that is independent of n. 

Thus, there exists u G L°°(0, T; H^) such that as n — t- oo (after the extraction of a subsequence) 



n 



XX 



u L°^{0,T;H^) weak*, 
"^u^^ L°°{0,T;H-^) weak*, 
^ Vu L°°{0,T;L^) weak * . 



Furthermore, since /i" — )• w in the sense of distributions (with v given by (2.5)), it follows from 
(2.9) that /i" r ^ vr in L°°{0,T; H^^) weak *. The equation for m" then yields u" ^ Ut in 
L°°{0,T; H~^) weak * and we conclude that u satisfies the first equation in (2.2). Uniqueness 
is a consequence of the linearity of the equation. It is possible to deduce additional regularity. 
Indeed, since u e L'^{0,T;H^) and ut G L'^{0,T; H^^), then , (cf. [23l Chap. 2]) we can assume 
u € C([0, T]; L^), which completes the proof of Theorem 



2.1 



3 Linearized stability for the Schrodinger— Burgers system 
3.1 Stability statement 

We now turn to the linearization of the full system ( |1.2[ ) (with e > 0) and consider the reference 
solution (e'''*r£, (^g) given in SectionjlJ cf. (1.8), (1.9). After a routine calculation, we arrive at the 
following linearized Schrddinger-Burgers system 

iut + Uxx ^ [ipe + b - 2£ rl)u - er^u + vre, 
vt + 2{ip^v)a; = 2e Re{rsu)x, 

supplementd with the initial data 

u{x,0) = uo{x) e iJ^(M), v{x,0) = vo{x) e i/^(M). 

We introduce the following notion of weak solution. 



(3.1) 



(3.2) 



Definition 3.1. Consider the reference functions if^ and in (1.8)-(1.9) for some e > 0. A pair 
(u, v) e L°°{0,T; x H-^(m.)) is called a weak solution to the Cauchy problem ([3l])-(|3^ if the 
following conditions hold: 

1. M is a weak solution to the Schrodinger equation 



lUt 



{ipe + b - 2er^)u - er^u + 



7 



2. v€ L°°{0,T;H^^ 



has the form 

v{t,x) ^ v{t,x) + ^'(t)(5s, 



(3.3) 



where v € L°°(0, T; ^^(k)), E := {{t^x) : x = 0,t G (0,T)}, is the Dirac measure on E, 
'i'(t) € L°°(0,r), and in the sense of distributions on M x (0,T) 



(3.4) 



This definition agrees with the one in Section |2.2| and, as far as the hyperbohc equation is 
concerned, generahzes the one in |20| for hyperbohc equations (but with a non-vanishing right- 
hand side). 



vt + 2{(p^v)^ = 2eRe(reu) 

with ipeSs, '■— 0. 
3. FinaUy, the initial data are assumed in a natural sense. 



Remark 3.2. Of course, under the assumption that v has the form (3.3) one can directly rewrite 
the equation for v as vt + 2{ipgv)x — 2£ Re(reit)j;, so that the condition ^Pe^T. '-^ need not be 
stated explicitly. Furthermore, from our proof of Theorem |3.3| below it follows that a solution v 



satisfying Definition 3.1 is actually the unique duality solution (in the sense of [S]) to the second 
equation of (3.11 (as follows easily from Theorem 4.3.2 therein, by taking the source-term into 
account). 

We are now in a position to state the main result of this paper. 

Theorem 3.3 (Linearized stability property for the nonlinear Schrodinger-inviscid Burgers sys- 
tem). Consider the reference functions r^ = ri.{x) and = ^Pe{x) given in Section^ cf. (1.6)- 



(1.7). Then, the solution {e r^,ip^) to system (1.2) is linearly stable in the sense that, for any 



T > 0, and any initial data (2.3), the system (3.1) admits a unique weak solution {u,v) = 
{u{t, x),v{t, x)) in the sense of Definition 3.1 with 



where Gt 



ll(u,w)llL~(0,T;_ffix//-i) < GT(||(u(0),u(0))||i/2x_ffi 

is a continuous function vanishing at the origin. 



(3.5) 



3.2 Proof of stability 

This section contains a proof of Theorem |3.3[ that is, a proof of well posedness for the problem 



(3.1)-(3.2|. To this end, it is necessary to carefully study the coupling between the two equations 
and deal with the lack of regularity of the coefficient of the second equation. We can first establish 
a well-posedness result for a regularized problem with smooth coefficients, as now stated. 

Proposition 3.4. Let Oj — aj{x) (j = 1, . . . ,5) be real functions in W'^'^{M.) such that a^ = 03 



up to a multiplicative constant, and fix any data uq G H'^{M.) and vq G H^{I 

iut + Uxx — CLiu + a2U + a^v, 
Vt + {aiv)x = (05 Reu)^,, 



Then, the problem 
(3.6) 



with prescribed initial data (w(0),f(0)) = {uq,Vo), admits a unique solution satisfying for allT > 
{u,v)eC{[0,T];H^)nC^[0,T];L^)xC{[0,T];H^)nC\[0,T];L^). 



Proof. In order to construct the local-in-tinie (strong) solution to (3.6), we follow the technique in 



24J[T5] and introduce an auxiliary system with non-local source, which can be tackled via Kato's 



theory [19]. This is necessary in order to write the system (3.6 1 without derivative loss (see [l5] 
for details). In other words, by refering to [15j for the motivation, we consider the linear system 



iFt + Fxx = aiF + - asOiVx 
Vt + a^Vx + a'^v = (05 Re u)x, 



a^a'^v - 03(05 Re u)j;, 



(3.7) 
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where F is the complex conjugate of F and F, u are defined by 



/ F{x, s) ds, 
Jo 



u{t, x) — Uq(x) 

u{t, x) ^ {A — 1)^"'" ((fli — l)u + a2U + a^v — iF) 
while the initial data read 

-2 



F{;0) = Fo e L\R), 
Once we have a solution 

F e C{[0,T]; L"^) n C\[Q,T]; H-^) 



v{-,0)=vo G H^R). 
veC{[0,T];H^)nC{[0,T];L^) 



(3.8) 

(3.9) 
(3.10) 



to the problem (3.7), one can argue (as in [121 Lemma 2.1]) and show that {u,v) given by (3.7)- 
(3.8) is actually the desired solution in the statement of the proposition. 

At this stage, we only sketch the argument since it is quite similar to the one in [HITS]. First, 
we write (3.7) as a system of three equations, by decomposing F into its real and imaginary parts, 
which allows us to obtain the abstract form 

Ut + A{U)U ^g{t,U), C/(-,0) = [/n, (3.11) 

with U ~ (Re F,lmF,v) and some initial data Uq. The key point is to decompose the operator 





A{U) = 



A 
-AO 

a^dx -\ 



in the form SA{U)S ^ = A{U) + B{U) for some operator B (again, see [1] [15] for a similar 
statement). In the present setting, we obtain such a decomposition by setting 



1 - A 
1 - A 




(1- A)i/2 



Note that S' : F X is an isomorphism, provided Y := {L'^f x and X (i?-^)^ x L"^ . The 
relevant properties satisfied by S (in particular about (1 — A)^/^) can be found in Section 8]. 
Observe that the right-hand side of (3.11) is linear in [/, so that it is straighforward to derive the 
necessary estimates for the source g and we may finally apply [191 Theorem 6] and conclude with 
the existence of F,v satisfying (3.10). 

We have thus established a local existence result for the Cauchy problem (3.6) and, now, in 
order to show that the solution exists for all times T > 0, we must derive certain a priori estimates. 
We begin with an energy inequality satisfied by solutions to (3.6). Multiply the first equation in 
(3.6) by ut, take the real part, and integrate over K, so that 



d 
dt 



dt 



dt 



dx -\ — - / ai|u| + Re / c?x = — 2Re / a^vutdx. 



(3.12) 



Since by assumption 
find 

Re / a^vut dx - 



ca^ for some constant c, we can use the second equation in (3.6) and 
a^vtu dx 

(a3u)2,a4U dx 

— / a^v'^ dx. 
dt .hs 



Re 


d 
di 


/ a^vudx — 
Is. 


Re 


Re 


d 
di 


I a^vu dx + 


Re 


Re 


d 
dt 


/ a^vudx — 

m 


Re 


Re 


d 
di 


1 a^vu dx — 
m 


Re 



; / 03 Reu(a3 Reu)a; 
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Plugging this result into (3.12), we arrive at the energy conservation property 



d 
dt 



+ ai\u\ +Rea2M + 2'Bji a^vu — a^v \ dx = Q 



and an immediate consequence is thus 



WUxWl < C + c||w||2 + C||u| 



(3.13) 



(3.14) 



Now, multiplying the second equation in (3.6) by v, integrating over M, performing various 



calculations, and finally applying Gronwall's lemma (assuming v smooth without loss of gen- 
erality, since the corresponding norm does not appear in the final estimate), we find ||w||2 < 



c + c{t) Jq \\ux\\2 + \\u\\2ds. Along with the estimate (3.14) and again by applying Gronwall's 
lemma, this gives 

Ml < c{t) + c{t) f\\u\\lds. (3.15) 
Jo 



Next, multiplying the first equation in ( |3.6[ ) by u, taking the imaginary part, and integrating in 
M, we obtain 

-— / \u\^dx = —lni / a2{u)^dx + lm / a^vudx. (3.16) 

2 dt Ib Ib Im 



Thus, an integration over (0,i), Gronwall's lemma, and (3.15) yield us 



t ft 

2||2 ,„ ^ „u\ , „U\ / Il„,2|l2 



\\ur2<c{t)+c{t) \y\\ids<c{t)+c{t) Wu'Wids 

Jo Jo 



and, with Gronwall's lemma again, < c{t). From (3.15), we see that ||w||2 < c(t) and thus, in 



combination with (3.14), we find 



(3.17) 



\K\\l + \Ml + \\v\\l<c{t) 

for some continuous function c = c{t) . 

To complete the proof, we need to an estimate on Vx,vt, Ut and Uxx and we only sketch this 
direct argument. First, we differentiate the second equation in (3.6) with respect to x, we multiply 
by Vx and integrate over M. We obtain ||wx||2 ^ c||wa;x|l2; together with terms that have already 
been estimated. Considering the first equation in (3.6), we multiply by Uxx, integrate over M, and 
use Young's inequality; this leads us to control Hi/JTilf ^ cHuflU -l-c. Now, by multiplying the first 
equation in (3.6) by ut, taking the imaginary part, integrating over M, and using Gronwall's lemma 

]+c for some small 7. These estimates combined 



< 



l\\Ux 



and Young's inequality, we control 
together give us precisely the regularity stated in the proposition which is now completed. 



□ 



To be able to apply Proposition |3.4[ we now introduce a regularized version of the linearized 
system (|3.1|. For S > 0, we define (p^ from ip^ by convolutio n of the sign function with a standard 



moUifier. (In particular, (p^{0) = 0.) Using Proposition 3.4 we obtain a solution {u' 



C([0, T];H^)nC\[0, T];L^)xC{[0, T];H^)nC\[0, TY.L'^) to the Cauchy problem (e being omited) 



iut + Uxx = + b — 2er^) u - 
Vt + 2{Lp^v)x = 2e Re{ru)x, 



er u 



(3.18) 



with initial data 



u\x,0) = ufi{x) e H^{R), v\x,0) ^vo{x) e H\ (3.19) 
With this solution in hand, we need estimates that are uniform in 6 and, clearly, the estimate in 



the proof of Proposition 



3.4 



ip^ approaches a discontinuous function as i5 
limit 5 — > 0. 



are not suitable, since they involve norms of derivatives of ip^ . Since 
0, these estimates provide no information on the 
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Lemma 3.5. For S > 0, th e solu tions {u^v^) E C{[0,T];H^) D C^{[0,T]; L^) x C{[0,T];H^) D 
C^{[0,T];L'^) to the system (plS]) satisfy 



e L°°{0,T;H^), uf e L°°{0,T-H-^), 
v'^ e L°°{0,T;H-^), vf £ L°°{0,T;H-^), 
ui^£L°°{Q,T-H'^), i;V* e L°°(0,r;L2) 



(3.20) 



uniformly in S > 0. 



Proof. In what follows, the constants c may depend on t. From (3.16), by using the duality 
l(/,5)l<ll/IU-i|l5lUi,wefind 



2<c + c/ ||u||2(ii + c / jif-i (it 
Jo Jo 



and, from the energy conservation property (3.13), 



Ikxlla <c + c\\uf2+c\\v\\H-i\\u\\„i+c I \(p^\v^dx 



(3.21) 



(3.22) 



Now, by multiplying the second equation in (3.18) by a test function in and integrating in 
time, we obtain (since C H^^) 



v\\h-i < c + c 



< c + c 



\lP^\v'^ dxdt + c I \\ux\\2dt 



(3.23) 



The latter inequality is a consequence of {<f^)^ < c\^p^\, which follows directly from (1.8)-(1.9). 
Also, the same equation yields immediately 



\\dtv\\'ij-i < c I \(p^\v^dx + c\\u,\\l 



(3.24) 



In view of the estimates (3.21 )-(3.24), it is necessary to bound Jj^ dx, only. Namely, by 

multiplying the second equation in (3.18) by w|</3''|, we obtain 



1 d 
2dt 



v^\f'^\dx + 2 I {vip'')^v\ip''\dx ^ I 2e(Rc ru)^v\ip''\dx 



and, since sgn(p^{x) = — sgnx (thanks to our choice of regularization). 



{v^%v\ip^\dx^ - / {vip^)^v\ip^\dx 



Vip\v\ip^\)or: dx 



1 



{vLp )^v(p dx - - I {vtp )xVtp dx 



vifi^ {vLp^)x dx + - I vip" {vip'')x dx = 



Thus, using again {ip ) < c|(p |, we arrive at 



v'^\ip^\dx < c + c I \\ux\\2dt + c 



(3.25) 



The estimates (3.21 1-( 3.25) and Gronwall's lemma lead to the desired properties in (3.20). The 



estimates for Uxx and Ut in ( 3.20 ) follow from the equation satisfied by u and the proof of Lemma 3.5 
is now completed. □ 
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Returning to the proof of Theorem 3.3 we see that (3.20) imply that there exist u,v,£, such 
that (for a subsequence) 





u 


in 


L 




V 


in 


L 






in 


L 



\0,T;H-^), (3.26) 
The distribution v obtained in this way is now checked to be a solution to the second equation in 



(3.1). As observed in Section |2.2[ this is not a trivial question since v has a Dirac singularity on 
the line S = {{0,t),t G (0,r)}, and (p is discontinuous on E. Following Section 2.2 it is necessary 
to use the information provided by the equation to define a suitable notion for the flux ipv. 
Our first observation is as follows: if is an open set contained in / x (0, T) with / n {0} 



then v\n in fact belongs to L°°{0,T; L^(r2)) and, moreover, v satisfies the second equation in (3.1 1 
in n. To check this, note that, in O, the function tp is smooth and bounded away from zero. 
Therefore, from the last property in (3.261, we find C/f in L°°{0,T; L'^^fl)) and, along with 
the second convergence property in (3.261, we obtain v\fi G L°°(0,T; L^(r2)) and ^ = v(p almost 
everywhere in 17. Since, according to (3.26), ^ satisfies the equation vt + 2^^: = 2e(rReM)a;, we see 
that 

vt + 2{ipv)^ ^2e{rReu)^ (3.27) 



in the sense of distributions in J7. In other words, the equation (3.27) is satisfied outside the 
singular line E. 

Given any point {t, x) with x ^ 0, one can see using elementary techniques (integrating along 
the characteristics) that the equation (3.271 can be uniquely solved with initial data vq & if^(K), 
since no characteristic ever crosses the line E and the right hand side is an function. This 
provides us with a function v defined everywhere except on E, and such that the one-sided limits 
v{0+,t) and v{0—,t) exist for all t > 0. This shows that v and v coincide outside of the singular 
line E. 

Let us now determine the equation satisfied by v in the whole of (0, T) x M. Let be a test 
function defined on K x (0,T). By computing (vt + 2{ipv)x,(l)) in the sense of distributions and 
after some easy calculations, we find 



where 



Vt + 2{ipv),x = 2eirReu)x + J{t)S^, 



J{t) = 2(p(0+)w(0+,t) - 2(p(0-)w(0-,t). 



Suppose that v has the form 



(3.28) 
(3.29) 

(3.30) 



that is, f — w has the form above and is not a more general distribution. 

According to the discussion in Section |2.2| and provided we define the distribution product 



(fSs := 0, 



(3.31) 



we find (in the sense of distributions) 



Vt + 2{(pv)x ^vt + 2{ipv)a; + ^'{t)Ss + 2ip^{t)5j: 
= 2e{rReu)x + {J{t) + ^''(t))(5s. 



So, V satisfies the equation (3.4) under the condition (3.31) if we set ^'(t) — — J*J{s)ds. Thus, 
we have established that ^ 

v = v~ I J{s)ds, (3.32) 
Jo 



where J is given by (3.29). It only remains to show that (3.30) does hold. As seen, v coincides 



with V outside the line E, and so the delicate point is establishing that there exists \E' G L°°{0,T) 
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satisfying (3.30). Namely, it is not clear that v\y: is a bounded function of i G (0,T), and not a 
more general distribution. (One could conceivably have, for instance, w|s = '^t=to for some tg.) 
The following technical lemma is in order. 

Lemma 3.6. Let w G L°°{0,T; H-^{R)) be given such that suppw C S = {(0,t),t G (0,r)}. 
Then, there exists w G L°°{0,T) such that w — w{t)ds- 

Proof. We need to find a bounded and measurable function w on (0, T) such that 

{w,(l){t,x)) = f w{t)(t){0,t)dt={w(S)Sx=o,4>{t,x)) 
Jo 

for every test function <j) E !D((0, T) x M) (the space of all compactly supported smooth functions). 
First of all, w is a distribution of order at most 1, as this follows from 



\{w,^)\< / \{wit),(f>)H-i^mdt 



< 



|w||L~(O,T;H-i)ll0(-,i)l|Hi(R) '^^ - C||'/>||ci((o,T)> 



with C depending on the support of (f), only. Therefore, according to Hormander [171 Theo- 
rem 2.3.5], the distribution w has the form 



(w,(?!))d'((o,T)xR) = (w^O>x(0,i))D'(o,T) + ('i«l,0(O,t))B'(o,T), 



(3.33) 



where wq G T>'{0, T) is of order and wi G ©'(O, T) is of order at most 1. We will show that, in 
fact, = and that wi is actually the bounded function w of interest. 

We deal first with the first term in the right-hand side of ( 3.33 1 and fix any 0(t, x) G D((0, T) x 
E) such that (t>{Q,t) = for ah t G (0,r). In that case, (|3.33) becomes 



\{w,4>) 



D'((0,T)xI 



4>.Mt) dwoit) 



(3.34) 



since a distribution of order is a measure. On the other hand, from the assumption w G 
L°°(0,r;F-i), we have 

|(u;,(/))25'((o,T)xR)l < c||(?;'||Li(0,T;_f/i)- (3.35) 



We will now construct a family of test functions for which ( 3.34 1 and (3.351 lead to a contradiction, 
unless Wq = 0. 

Define a family of functions {On{x)}n>o in H^{M.) by 



-l~x, x G (-1,-1/2), 

X, x G (-1/2,1/2), 

1-x, XG (1/2,1), 

otherwise, 



and set On{x) := ^6Q{nx). For all n, we have 



,(0)=0, 0^(0) = 1, \\9n\\HHm< 



(3.36) 



Now let (t)n{t,x) — 9n{x)'il>{t), where ip is a. test function on (0,T). From (3.34) and (3.36), we get 



\{w,(f'n)\ 



i'it) dwoit) 



(3.37) 
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On the other hand, from (3.35 )-(3.36 1 we find 



\{w,K)\ < / \m\-dt. 



(3.38) 



The estimates (3.37) and (3.38) give /g 'ip{t)dwo{t) = 0. Since ip was arbitrary, we conclude that 
Wo = 0. 



We now claim that the distribution wi defined on (0, T) in ( 3.33 ) is, in fact, a bounded function. 
Let Ip he a test function on (0, T) and let (f> = 4'{x) be a test function on M such that 0(0) — 1. 
Since wq ~ 0, by using (3.33) we find 



\{wi,^)\^\{w,cj,{x)^{t))\= / {w{t),cl){x))H-ixmmdt 

Jo 

< \\'w\\L-=-{0,T;H-^)\\(f>\\HA\M\L^a,T) = C"]] ^|| (o,T) ■ 

By a density argument, we may suppose that the previous estimate is valid for any ip e L^{0,T). 
Therefore, wi is a continuous functional on L^(0,T) and so is in L°°{0,T). The conclusion of the 
lemma is precisely (3.33) with w = wi. □ 



To complete the proof of Theorem [33] we need to establish uniqueness. Consider two solutions 



{ui,vi) and {u2,V2) satisfying the assumptions in Definition 3.1 Denote u and v the differences 
ui — U2 and vi — V2- Since the system (3.1) is linear, u and v satisfy (3.1) and so uniqueness 



of solution will follow if we show that any solution with vanishing initial data does vanish on 
(0, T) X M. For this, the crucial estimate is 



1 



dt 



^ da; < 2e 



Re{ru)x\\(pv\. 



(3.39) 



From now, we omit the subscript e. Suppose that w(0) = and note that w(0) — f (0). With the 



estimate (3.39) in hand and from Gronwall's lemma and the fact that u € L°° {0,T; H^), we find 
v{t,x) = 0. In that case, (3.29) and (3.32 1 imply that, also, v{t,x) = 0. Then, the first equation 
becomes iut + u^x — {f + b— 2er^)u — er'^u, with u G L°°{0, T; H^). This allows us to multiply it 
by u and obtain the conservation of mass (3.16) without the latter term. Gronwall's lemma and 
m(0) = immediately give us u i 



0. 



It only remains to check the estimate (3.39). From the equation (3.28) and after multiplication 
by ipv and integration over (— oo, 0), we obtain 



1 



lpipdx + ip^{0-)d^{0-,t) = 



2e Ke{ru)x^v dx 



and, since ip > on (— oo, 0), 



-dt / v^(pdx< / 2e\Ke{ru)xfv\ dx. 

^ J —oo J —oo 



Similarly, integrating over (0, +oo) and using (/3 < 0, we find 



zdt 



+ 00 /' + C30 

rr2 



V (fidx < 



2e\ Iie{ru)x^v\ dx. 



Summing the two previous estimates gives (3.39). (Note in passing that (3.39), from which 
uniqueness follows, is satisfied since ip as an entropy-satisfying shock and would not hold for 
a "rarefaction-shock"; see ^I] for stability statements about this issue for general solutions to 
linear hyperbolic equations.) This establishes the uniqueness property and completes the proof of 
Theorem [331 □ 
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4 Numerical experiments 



4.1 Proposed setup 

In this section, we present some numerical computations which provide an illustratation of our 



theoretical results. We compute numerically the solutions of, both, the full nonlinear problem ( 1.2 1 



and the linearized problem (3.1). With these simulations in hand, we can check numerically the 
stability property (1.4 1. Recall that the linearization approach is based on the formal substitution 
V 1-^ Lp + 6vs and u e*''*r + Se^^^us, where {u,v) are the solutions to the full problem (1.2). 



This leads us to a linear system satisfied by the perturbation {us,vs) with initial data {us,vs). 



The claim is that the solution (u, v) to the full system ( 1.2 ) with perturbed initial data (urof (0) + 
Sus, Wrcf (0) + Svs) remains close to the perturbed exact solution {urcf + Se^^^ug, Vrd + Svg). 

Consequently, in order to numerically demonstrate the validity of the linearization procedure, 
we will proceed as follows: 



1. Compute the solution {ug^vs) to the linearized system (3.1) with some initial data (us,vs) 



2. Compute the solution to the full system (1.2) with perturbed initial data (ui-ct(O) + 

(5M5,Wrof(0) + 5V5). 

3. Compare (w, v) with the perturbed exact solution (wref + Se'^^'^ug, Vrei + Svs), and check that 
the difference remains bounded for any finite time. 

We will rely on a finite volume method and will not seek for an in-depth numerical analysis of 
the solutions, but rather visualize the linearized stability property of interest around our particular 
reference solution. 



4.2 Numerical results 



Specifically, the numerical approximations to (1.2 1 is based on the finite volume scheme described 
in [3^ . The Schrodinger equation is solved using a finite difference semi- implicit Cranck-Nickolson 



scheme, and a Newton algorithm allows us to deal with the nonlinear term 



Burgers equation 



in ( 1.2 ) is solved by a semi-implicit Lax-Friedrichs scheme. Concerning the linearized system (3.1 ) 



we employ the Lax-Friedrichs scheme for the second equation and a standard Crank-Nickolson 
scheme for the Schrodinger equation. 

Another issue is the numerical computation of the reference solution 7'g 



in (1 



We have 

employed an explicit Euler scheme after transforming the second-order equation (|1.8| into a system 



of two first-order equations. Since the Euler method is not particularly accurate, we rely on a fine 



mesh (with 20000 points) so that the solution to (1.8) provides a suitable approximation for our 
purpose. This is confirmed by testing the case e = 0, in which the solution (1.6 1,(1.7 1 is available 
explicitly. In our numerical experiments, we choose the parameter values e — 0.1,6 — —1.5 in 
(1.2), (1.8), and (1.9), and S — 0.1 in (1.3 1. The tests are performed with 20000 spatial points on 

0.002) and time increment r = 0.0001. 



For the 



the interval (—22, 22) (giving a spatial step h 
perturbations {us,vs) we choose a gaussian curve . 

In Figures [l] and [2] we present the solution {u,v) to the full system (1.2) with perturbed 
initial data (wrof(O) -I- 5us,Vrof(0) + 6vs) (dashed line), along with the perturbed exact solution 
(uiof +<5e'''*U5, Vycf+Svg) (full line) at the final time T = 1. Recall that the difference between these 
two quantities should not vanish, but only remain bounded, and it is indeed what we observe. 

In Figures [3] and [4] we display the solutions to the linearized system (3.1). Recall that the 
solution to the linear equation with discontinuous coefficient (3.1) contains a Dirac delta along 



a; = 0. Since our scheme contains some numerical viscosity, this Dirac delta is smoothed out and 
appears as a steep spike in the solution, which can also be observed in Figure [2j 
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Figure 1: Perturbed exact solution and solution with perturbed initial data - Schrodinger equation 



Comparison of perturbed exact solution and solution with perturbed initial data, Burger's equation 



Exact solution with perturbation 



Computed solution with perturbed initial data — • 
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-0.5 



0.5 



1.5 



Figure 2: Perturbed exact solution and solution with perturbed initial data - Burgers equation 
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Solution of the linearized problem, Schroedinger equation 
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Figure 3: Solution to the linearized problem - Schrodinger equation 



Solution of the linearized problem, Burger's equation 
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Figure 4: Solution to the linearized problem - Burgers equation. 
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